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Abstract 

RGEs for coefficients of dim-5 operators giving rise to neutrino 
masses in tlie seesaw mechanism are written down in tlie SM, 2HDM 
and MSSM, and solved numerically. RG evolution of these coefficients 
modifies tree-level seesaw predictions for neutrino masses and mixing 
angles in SO(10)-type GUT models as strongly as quark Yukawa cou- 
pling evolution. 
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The MSW solution of the solar neutrino problem |1|, Q has revived in the 
last few years interest in GUT model predictions for the superlight neutrino 
masses generated via the seesaw mechanism [§, ||. Such models predict the 
neutrino mass matrix (sometimes up to a normalization factor) using defi- 
nite relations between quark Yukawa couplings, Dirac type neutrino Yukawa 
couplings and Majorana mass matrix for the right-handed neutrinos. In this 
letter we elaborate on one technical aspect of such analyses: In comparing 
the GUT relations with low energy charged fermion and neutrino masses one 
has to perform a RG evolution of relevant parameters. The RG evolution of 
charged fermion Yukawa couplings from the weak (Myy) to the unification 
scale (Mx) is standard Q and has been included before [@, |^. The RG evo- 
lution of neutrino masses has so far been neglected or treated inconsistently 
[0]. We furnish this lacking element and demonstrate its possible numerical 
significance. 

GUT models we are interested in contain superheavy SU{2) x U{1) - 
singlet Majorana fields with a Dirac-type Yukawa coupling to the SM neutri- 
nos and Higgs boson. To describe physics far below the GUT scale we should 
decouple all heavy states and use an effective theory containing light fields 
only. Tree-graph exchange of superheavy right-handed neutrinos gives rise, 
after integrating them out, to dimension-5 operators coupling two left-handed 
lepton doublets and two Higgses. Couplings of this kind, being suppressed 
by the inverse heavy mass factor, are neglected in the standard approach in 
which one retains in the effective lagrangian renormalizable interactions only. 
In order to obtain the seesaw neutrino masses of the order 0{1/Mx) one 
must keep in the low energy lagrangian also those nonrenormalizable higher 
dimension terms. Similar couplings arise also in some string models and 
our RG analysis is also applicable in this case starting from the scale where 
the effective theory is SUc{3) x SUl{2) x f/y(l). We write down one-loop 
RGEs for the coefficients of the relevant dimension-5 operators which de- 
scribe their evolution down to M]y scale. When Higgs fields are replaced 
by their VEVs this coupling yields a superlight Majorana mass term for the 
left-handed neutrinos. In the Standard Model (SM) as well as the 2 Higgs 
Doublet Model (2HDM) we need only one dimension-5 operator - Oi . In 
the SUSY case we have to include other operators related to Oi by su- 
persymmetry transformations because they mix under renormalization. We 
solve numerically those RGEs in the SUSY case and find that they introduce 
modifications to the tree-level seesaw formula of the same order as the quark 
Yukawa coupling evolution. 

While evolution of the neutrino masses themselves does not provide very 
valuable information because the overall scale of heavy Majorana mass is not 
fixed precisely in GUTs, evolution of the neutrino mixing angles (e.g. if a 
definite texture for the Majorana mass matrix is assumed) may be important 
in comparing GUT predictions with MSW solar neutrino problem solution. 
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Superlight neutrino masses arise via the seesaw mechanism after diago- 
nahzation of the mass term (M ^ m, generation indices neglected): 



( V is the left-handed neutrino and n is the left-handed antineutrino). 

It has two eigenvalues: mh^avy ~ M and mught ~ m? /M, with the light 

neutrino: i^ught = cos a v + sin a n where sin a m/M. 

This viewpoint does not allow, however, any simple analysis of the neu- 
trino seesaw mass running. Moreover, a serious problem arise, in theories 
(like MSSM) in which electroweak symmetry breaking proceeds radiatively 
[^]. On one hand, the tree- level VEV of the Higgs boson(s) vanishes at Mx 
scale and no neutrino mass is generated there. On the other hand, to get 
Higgs VEVs one must use RGE of the effective theory with all superheavy 
particles decoupled. As a result neutrinos stay, at first sight, massless. The 
solution of the two above mentioned problems emerges when one retains in 
the effective lagrangian higher dimension operators generated by the decou- 
pling procedure. The heavy Majorana neutrino n with the Yukawa coupling 
in the GUT lagrangian of the form: 

AL = -Y^^n\i/iHj (2) 

( eij - antisymmetric SU(2) tensor; a,b - generation indices and Hi is the 
Higgs doublet transforming as (2, 1/2) under the electroweak gauge group 
which in the 2HDM and MSSM gives masses to the up-type quarks) when 
integrated out, gives rise to the nonrenormalizable term in the effective la- 
grangian 

A^e// = \ of Of (3) 

where 

of = {eikHil'l){ejiHjl\) 
and cf is a symmetric matrix of coefficients: 

cf{Mx) = 2 Y'^^iM^y^Y''^ 

(M is the Majorana mass matrix of n). It is easy to see that when Higgs field 
is replaced by its VEV, (3) gives rise to the neutrino mass term which up to 
0{1/Mx) is equal to the mass term given by (1). 

RGE for coefficient ci gets contributions from the proper vertex correc- 
tions due to gauge boson exchanges and Higgs self coupling, and from exter- 
nal wave function renormalization which introduces also the dependence on 
Yukawa couplings. In the 2HDM we get: 
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(4) 



( t = log{Q / Mw) and summation over the generation index c is 

understood ). 

In the SM there is an additional contribution to the RHS (because the 
Higgs boson couples to all fermions): 



6 tr (r,rj")+2tr (f^F/ 



„ab 



Where Y^, Y^, Yi are up-quark, down-quark and lepton Yukawa matrices 
respectively and A2 is the Higgs self coupling: 



(5) 



In the MSSM we have to introduce other operators related to Oi by 
SUSY because they mix under renormalization due to diagrams in which 
gauge bosons are replaced by the corresponding gauginos. The relevant part 
of the effective lagrangian reads: 
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AT ^O'b r^ab , „ab r^ab , „ab r^ab , „ab r^ab 

^^eff - ^ Ci Ui + C21 U21 + C22 U22 + ^ C3 (J3 



(6) 



with 



O2I = {eikhiLl){ejiHjli) 



O22 — i^ikhilDiejiHjLi) 



or ^ {e,,h,Ll){ejihjLt) 

/i is a higgsino and is a slepton. Operators Oi to O3 form a basis 
in the space of operators which mix under renormalization. 

For a superpotential containing the singlet neutrino superfield TV": 



w ^ - M^'^N^'N^ + Y^-^N^'eijLlHj 
we get after integrating out and its supcrpartner: 

Icf = C^J = 0^2 = = ycai^j^-X^cdydb 



(7) 



(8) 
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The RGE for cf' in this case reads: 



+ 
+ 
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(2 GmGii + 6 G2HG2L) {4i + C21 
(2 GihGil + 2 G2HG2L) (c22 + C22 
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In equation (9) Gil, G2L, Gih, G2H are U(l) and SU(2) gaugino 
couphngs to leptons and Higgs boson. Yim and YihE are Yukawa couphngs 
of leptons to Higgses and higgsinos respectively. In the strict SUSY limit they 
are all equal to the corresponding gauge and Yukawa couplings gi, g2 and 
Yi. In the same limit A2 = 91+92- We display, however, this equation in its 
more general form which allows one to decouple (when necessary) all heavy 
sfermions in one collective threshold. The RGEs for other MSSM couplings 
written in the same general form can be found in [^. 
The remaining RGEs read: 



±ab 
df' 
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In (10-12) we implement all SUSY relations between couplings but we re- 
tain the distinction between cf^ and Cg" because in principle one could have 
a model where they are not symmetric. In (4) and (9-12) there is no mixing 
with operators involving the other Higgs doublet (2,-1/2) because there 
is no appropriate coupling in the MSSM and simple 2HDMs. Contributions 
from Yukawa couplings come only from wave function renormalization and 
therefore do not cause operator mixing. This fact simplifies considerably the 
analysis of lepton mixing angle evolution presented below. 

In order to demonstrate the potential importance of the RGE evolution 
for neutrino masses and mixing anles we have solved equations (9-12) for the 
SUSY case numerically, evolving simultaneously gauge couplings and Yukawa 
couplings in one-loop approximation and including nonlinearities from the 
3'''^ family. We take a^iUz) = 0.12 and sin^ Ow = 0.233 and allow 
Ytop{Mx) to vary from 0.2 to 6. Our unification scale is Mx = 10^^ GeV. 
We consider two distinctly different cases: 
i) SO(10)-type unification for Yukawa couplings with 

Yt{Mx) = Y,{Mx) = Y^{Mx) 
a) An example of SU(5)-type Yukawa coupling unification with 

YtiMx) = 10 YhiMx) = 10 K(Mx) . 
After evolving the Ys down to Mz we fix tan/? fitting 171^-=!. 784 GeV: 
mT-{Mz) = {Yr{Mz) / 92) Mw cos P and calculate mtop{Mz) and mb{'mb) 
which is presented in Fig. la . Yt{Mx) as a function of rritop is presented in 
Fig. lb. This sets the background for the evolution of coefficients q and 
allows one to assess the proximity of the Landau pole. Our aim in this 
note is to demonstrate the effects of the neutrino masses and mixing angles 
evolution rather than performing a detailed study of Yukawa coupling unifi- 
cation. Therefore, we neglect possible deviations from the assumed relations 
between Yukawa couplings at Mx (threshold corrections) which could re- 
sult in mi, closer to its experimental value mi,{mi,) = 4.25 ± 0.10 GeV 
[|] [|. This simplified approach suggests (see ^ for a detailed analysis) that 
Yukawa coupling unification favours a heavy top quark. 

^We neglect also possible effects of decoupling superpartners at some scale Msusy = 
0(1 TeV) on the RG evolution. 
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We work in the approximation of small 3'''^ generation mixing. In evolv- 
ing Yt, Yb, Yt (or eg. Y^, Y^, Y^ ) we neglect this mixing altogether and in 
evolving quark and neutrino mixing angles we use approximate RGEs [|, § 
which for the MSSM read: 



(13) 



di 



Those simple RGEs arise because Yukawa couplings responsible for lepton 
mixing angle evolution appear exactly in the same way in (9-12). The mix- 
ings between different operators Oi - O3 and different generations can be 
factorized substituting c"^ = bi and, in fact, one could guess (13) with- 
out calculating (9-12) explicitly. The results are presented in Fig. 2a - 2d. In 
all figures there are 2 sets of curves: for Yukawa unification i) - solid lines 
and a) - dashed ones. 

Fig. 2a presents running of the coefficient ci for the 3'''^ and for the 1*^* 
(or 2"*^) generation. For heavier top, close to the fixed point of the RGEs, 
Yukawa contributions overcome the gauge coupling contributions and change 
the direction of evolution of coefficients q, just like in the case of quark 
Yukawa couplings. 

Fig. 2b presents running of the neutrino 1-3 (or 2-3) mixing angle, which is 
sizeable only for large Yr - i.e. in the case of SO(10)-type of Yukawa cou- 
phngs unification and large rritop- 

Fig. 2c presents the joint effect of both KM and neutrino mixing angle evo- 
lution. If we assume that they are equal at Mx as predicted by some 
models then their ratio at Mw, which is plotted in Fig. 2c, is due only 
to RG running . Running of the neutrino mixing angle is of the same order 
of magnitude as KM angle running and partly cancels effects of the latter in 
^is/^il^^- '^^^ that evolution even slightly strengthens the discrep- 

ancy between 6^^^ = (0.3 — 1.0) x 10"^ and the small angle MSW solar 
neutrino problem solution: 9'(^ = (3 — 5) x 10^^ [If. (In SUSY-GUT seesaw 
the solar neutrino problem is solved by i^r transition.) 

Finally in Fig. 2d we present the values of coefficient r] incorporating the 
effects of both neutrino and quark mass running in the seesaw formula: 



(no sum over generation index a) 

Tji for i = 1,2 are independent of rritop because in these quantities trace- 
type Yukawa coupling contributions to the evolution of Ytop and ci cancel each 
other. For rjs the trace-type contributions ( ~ triYtY^) ) to the evolution 
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of Ci and Yt cancel each other as well but the evolution of Yt due to the 
non-trace contributions ( ~ (YfYi) ) is stronger than analogous evolution of 
Ci ( ~ (Yt-Y^) ) and results in a large rj^ for a heavy top quark. Values of 
coefficients rj have been given in 0] (with light quark masses renormalized 
at lower scales). Even though numerical differences (~ 20%) are insignificant 
considering our ignorance of the overall scale of M we believe our results are 
more complete. We include carefully the effective dim-5 coupling running 
and state the assumptions concerning the charged fermion Yukawa coupling 
evolution more clearly, including assumed Yukawa coupling unification - the 
essential element of seesaw neutrino mass predictions. 

To summarize, superlight neutrino masses arising via the seesaw mecha- 
nism can be most naturally viewed as a manifestation of an effective dimension- 
5 operator. This approach is the only one available in models with radiative 
SU(2)xU(l) symmetry breaking. The RG analysis based on dimension-5 
operators provides a systematic and unambiguous method of incorporating 
radiative corrections to the GUT predictions for neutrino seesaw masses. 
Running of the coefficients q is an as important element of neutrino mass 
predictions from GUTs as quark Yukawa coupling running. For a heavy top 
quark the most important contribution to RGEs comes from Yukawa cou- 
plings of the 3'''^ generation. Any specific model predicting neutrino mixing 
angles from GUT seesaw must take into account equations (9-12) when it is 
compared with experiment. 

We thank Prof. S.Pokorski and Dr. M.Olechowski for discussions which 
initiated this study. Z.P. would like to thank Prof. D.Wyler for kind hospi- 
tality extended to him during his stay in Ziirich. 
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FIGURE CAPTIONS 



Figure 1. a) mb{mb) and b) Yt{Mx) shown as a function of rritop in 
unification schemes i) and ii). 

Figure 2. The result of the dim-5 coupling evolution: a) coefficients 
ci for the 1^* and 3"^^ generation and b) lepton mixing angle ^J'a- 
Joint effect of dim-5 coupling and quark Yukawa coupling evolution: c) 
e'(^{Mw)/0^-i^{Mw) and d) coefficient rj for the 1"* and 3'"'^ generation. 
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